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EXPERIMENTAL AND ANALYTICAL STUDIES OF A MODEL HELIOOPTER ROTOR IN HOVER 


F. X. Caradonna and C. Tung 
Aeromechanics Laboratory 

U.S. Army Research and Technology Laboratories (AVRADCOM) 

SUMMARY 


The present study is a benchmark test to aid the development of various rotor 
performance codes. The study Involves simultaneous blade pressure measurements and 
tip vortex surveys. Measurements were made for a wide range of tip Mach numbers 
including the transonic flow regime. The measured tip vortex strength and geometry 
permit effective blade loading predictions when used as input to a prescribed wake 
lifting surface code It is also shown that with proper inflow and boundary layer 
modeling, the supercritical flow regime may be accurately predicted. 


SYMBOLS 


A ratio of vortex circulation to maximum blade-bound circulation 

& aspect ratio 

C|^ sectional lift coefficient 

d radial distance from a vortex to a flow-field point 
R radius of the rotor blade 

r radial distance from the rotor center of rotation 
vortex-induced velocity 
Vj^ residual velocity in the wake 
y r/R, nondimensional radial coordinate 

z axial distance from rotor 

i; rotational speed 

Y azimuthal angle measured from the point of blade overhead passage 

Vy vortex age, the azimuth angle, V, when vortex strikes the probe 


^Presented at the Sixth European Rotorcraft and Powered Lift Aircraft Forum, 
September 16-19, 1980, Bristol, England. 



1 . INTRODUCTION 


The past two decades have seen a continuing develoiment of methods to predict 
rotor hover performance with increasing accuracy. These methods Include lifting line 
(refs. 1-3), lifting surface (refs. 4-6), and finite difference (ref. 7) methods. 
Practically speaking, none of these methods is self-contained; they all require the 
specification of empirically obtained wake data (strength and geometry) in order to 
have a correct downwash distribution. Inevitably, the development of these codes 
becomes a tuning process in which it is determined just how detailed and accurate a 
wake description must be. This stage of code development places great reliance on 
the available body of experimental rotor data. 

The available rotor data include a sizeable number of tests where detailed blade 
loading is obtained using surface pressure transducers (refs. 8-11), and more recently 
by laser doppler velocimetry (refs. 12-14). There is also a number of tests in which 
the rotor wake geometry is defined by flow visualization techniques (refs. 3 and 5) 
for a wide variety of blade configurations. Of the various wake studies, only 
Boatwright (ref. IS) and Cook (ref. 16) made detailed investigations of the wake flow 
structures. Cook's work is especially significant in that he was able to measure the 
strength of the tip vortex by a curve-fitting technique using hot-wire data. However, 
there seem to be no useable data in the literature in which simultaneous blade load 
distribution and wake measurements are made. 

It is the xntention of the present study to help fill this gap in the literature. 
This paper will describe the experimental set-up in which steady blade pressures were 
obtained using hub-mounted transducers and tip vortices were measured using Cook's 
technique. The data obtained are for unstalled flow ranging from the low subsonic to 
transonic conditions. It is shown herein that the measured wake geometry differs 
significantly from previously published low-aspect-ratio data (ref, 5). This differ- 
ence is reflected in an inability to correctly predict the measured blade loading 
(using Summa's prescribed wake lifting surface code (ref. 6)) when this classical wake 
geometry is used. 


2 . THE EXPERIMENT 


The data presented in this paper were gathered in the Army Aeromechanics Labora- 
tory's hover test facility, a large chamber with special ducting designed to eliminate 
room recirculation. The rotor, situated in the center of the chamber, was mounted on 
a tall column containing the drive shaft (fig. 1). The rotor employed two cantilever- 
mounted, manually adjustable blades with half degree precone. These blades used an 
NACA 0012 profile and were untwisted and untapered. An aspect ratio of 6 was chosen 
in order to maximize Reynolds Number and available instrumentation space. The blades 
were grooved to accommodate 60 pressure tubes each. These tubes connect to a special 
cluster of three 48S8 Scanlvalves (using Statham PA 856-15 transducers) driven by one 
SS5-48 solenoid drive mounted in the rotor hub. This arrangement permits an ample 
number of ports for five measurement locations — three radial locations on each blade, 
with one location being identical on both blades for comparison purposes. The 
Scanivalve stepper motor was actuated by a digital data system which acquired the 
data, computed the centrifugal pressure drops, and displayed the final pressure 
distribution. After manually adjusting the two blades, the pressure data was also 
used to check the equality of loadings. The pressure data at the 0.8 R radial station 
are compared for the two blades in figure 2. No significant differences in the 
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loadings were seen for any operating conditions. (Additional indication of this 
loading equality is that no consistent difference in the two shed vortices was 
found.) The resulting pressure distributions for collective pitch settings of 5*, 

8”, and 12“ are shown in figures 3, 4, and S. These and other pressure distributions 
are tabulated in appendix A. It is seen in figures 3, 4, and 5 that the inboard 
pressure distributions are only slightly affected by rotor speed. However, the 
outboard sections show conside: able pressure alteration and shock development as the 
tip Mach number approaches near sonic values. Overall, however, the spanwise load 
distribution (obtained by pressure integration) is remarkably little affected by tip 
Mach number (fig. 6). In addition, the tip pressures were compared with those of 
reference 11 and are seen in figure 7 to be nearly identical. 

Wake data were acquired with a traverse-mounted DISA 55P01 hot-wire probe mounted 
beneath the rotor. The probe was oriented with the wire being tangent to the rotor 
tip path. This permits measurement of the magnitude of the vortex induced velocity 
when the remainder of the rotor downwash is properly accounted for. It also 
excluded the effect of axial velocity on the induced velocity measurement. Dr.ta from 
the wire are acquired at various points along the tip vortex trajectories and can give 
both the tip vortex geometry and strength. One problem with this approach is that the 
vt ;tex trajectory is not steady and the probe location (which is chosen by an on-the- 
spot decision as to where the number of vortex core "hits" is maximized) contains some 
as yet undetermined error. The resulting data stream has considerable variability. 
However, in order to be certain of the vortex location, the only acceptable data are 
chose where the vortex core actually hits the prcbe. In the digitization process 
(done off-line at a reduced tape speed), the above-mentioned data system was coded to 
look for ana accept on’y those data which showed the characteristic signal dip 
(wherein the minimum velocity is very close to the vortex translation speed) which 
indicates a probe-vortex strike. This tarns out to be a very small percentage of the 
total amount of data actually recorded. A typical hot-wire trace displaying Che 
above-mentioned variability is shown in figure 8. 


3. HOT WIRE DATA ANALYSIS 


The idea of the current data analysis is that a tip vortex should look like an 
infinite line vortex to a sufficiently close probe. Unfortunatly , the probe mea- 
sures not only the velocity induced by the vortex at hand, Vj^, but also that induced 
by the blade and the remainder of the wake system as well, Vjj. The problem in analyz- 
ing the probe data is, then, how to separate this residual velocity, Vj^, from the 
immediate vortex-induced velocity, Vj. Cook (ref. 16) handled this problem by assum- 
ing that the residual velocity was constant and given by the translation velocity of 
the tip vortex. He then was able to find the vortex strength by a fitting process. 
This strength was found to be far less than the computed maximum blade bound circula- 
tion of the single, full-scale blade used in that test. It was also found that Che 
vortices measured were distinctly nonclassical in that they contained a large rota- 
tional region outside of the viscous core. In what follotra, we shall use a process 
very similar to Cook's in analyzing wake data. 

First consideration is given to the vortex trajectories. Figure 9 shows the 
•ixial and radial components of the vortex trajectories for a pitch setting of 8*. 

This figure gives data for a wide range of rotor speeds, and it is apparent that the 
trajectory is essentially independent of tip speed — even into the transonic regime. 
Figure 9 together with figure 6 suggests that the nonlinear transonic flow on Che 
blade has little effect on the far-field induced flow as long as the local lift is 
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not greatly altered. Also plotted on this figure Is the vortex trajectory given by 
Kocurek's wake-fitting formula for rotors in free air. Although the axial component 
of the trajectory compares well with Kocurek's formula, there appears to be a greater 
discrepancy in the contraction than can be explained by measurement error. The vortex 
trajectories for pitch settings ranging from 5” to 12" are given in figure 10. 


The present aim in analyzing the rotor wake is only to find the vortex strength 
and not a complete description of the structure. This strength will be found by 
fitting the wake data to the velocities obtained from an appropriate combination of 
invlscid, two-dimensional vortices. The velocity from one such vortex is given by 


fiR “ ^ 4 ttA d/R 


( 1 ) 


where the strength of the vortex is described by A, the ratio of the vortex circula- 
tion to the maximum bound circulation of the blade. (This could be determined by the 
pressure data because the circulation peak is not very sharp and is quite close to 
Che tap location.) To accomplish this fitting, it is first necessary to convert the 
spatially dependent equation (1) into a time-dependent expression, as the vor ex data 
are time-based. Assuming that A is constant (which seems to be true within reason- 
able error bounds) , the conversion to a time-dependent function is accomplished by 
expressing d as a function of time using the vortex trajectory data of figure 9. 

The next step is the determination of the residual velocity, Vr, which must be vec- 
torially added to before a comparison can be made with the probe data. We have 

done this in two different ways: 

1) The first way to determine involves very young vortices (about 50° old). 

For these it was assumed that Vj^ was given by the vortex trajectories (fig. 9). The 
fitting process always commenced when the vortex core hit the probe and ended when the 
following blade passed over; this assured the simplest possible flow field, as there 
would be vortices on only one side of the probe and minimal influence of vortex sheets 
and blade bound vorticlty. Figure 11 shows some typical con arlsons of probe data with 
the fitting expression. This figure shows the vortex velocity-time traces for pitch 
settings of 8° and 12°. It is seen here that the fitting curve provides a good match 
to the data outside of the immediate core region. Furthermore, the vortex strength 

is very close to the maximum blade bound circulation. 

2) A second means to determine was required in analyzing older vortices 

(about 210° old). The flow is more complex in this case, as the probe always lies 
between two vortices in the fitting region, and the expression for the vortex- induced 
velocity is correspondingly complicated. In fact, for this case was determined 
using three vortices — one outboard of the probe and two inboard. Again, the data 
were fit for the time period between a probe-vortex strike and the subsequent blade 
passage. It was found that with Vp determined by the vortex trajectory data, it was 
not possible to obtain a good fit of the classical vortex expression to tne wake data. 
Instead, we found that a better value for Vj^ was found by use of the minimum mea- 
sured velocity between two vortices. At this point, the vortex-induced velocity is 
i-mall, but not zero (due to the differing instantaneous translation velocities of the 
three vortices). The minimal Induced velocity is calculated (assuming some value of 
A) and subtracted fiom the minimum measured intervortex velocity to obtain Vr. This 
task was rendered quite simple by the fact that the radial component of these veloc- 
ities turns out to be very small (this was checked by calculations and measurements 
with a second probe) . Since the two methods above do not give the same value for the 
residual velocity, it is clear that is not a constant in this case. We assume. 
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however, that it changes sufficiently slowly to render the fitting process meaning- 
ful. In fact, the results thus obtained are consistent with the young vortex data. 
Figure 12 shows some typical comparisons of the older vortex data with the fitting 
expressions. This figure shows the fittings for pitch settings of 5**, 8”, and 12". 

It is seen that the 8" and 12" cases show vortex strengths which match the maximum 
blade-bound vorticity very well. At 5" pitch, ho%rever, the strength is seen to be 
considerably smaller. 

It seems from the above data, which are taken at a low rotor speed, that the tip 
vortex develops its full strength very early in life (mainly before 50"). Although 
there is a fair amount of variability between vortices, it is rather striking that 
very many vortices closely approach a classical Rankine vortex in appearance. Fur- 
thermore, the vortices (except for the 5" case) seem to contain all of the blade 
circulation. This vortex strength and structure differs markedly from the result 
obtained by Cook and probably reflects the considerable differences in blade geom- 
etries. As rotor tip speed Increases (fig. 13), however, there appears to be an 
Increasing departure from the Rankine vortex appearance. Nevertheless, the nondlmen- 
slonal vortex strength seems unaffected by tip speed. 


4. COMPARISON OF THEORY AND EXPERIMENT 


In order to Integrate the present wake and loading data into a believable whole. 
It Is necessary to be able to reproduce the blade loading computationally. He have 
chosen to do this using A.M.I.'s lifting surfac.: code (ref. 6). This is a very 
flexible, compressible, lifting surface code which can handle either prescribed or 
free wakes. 

Initial efforts to compute the blade loading were done using the Kocurek wake 
geometry (ref. 5). The resulting computed thrust coefficient was too high by about 
20%. The next step was to compute the loading using the measured vortex locations and 
strength. Figure 14 shows a comparison of the measured and computed loading using the 
measured vortex parameters for a collective pitch of 8" (the trajectory is given by 
fig. 9 and we choose A - 1.0). The comparison is now considerably improved and the 
thrust coefficient is overpredicted by less than 5Z. In view of the previously men- 
tioned uncertainties in the vortex trajectory measurements, these computations were 
also performed with the entire vortex trajectory perturbed such that at 4* ■ 180", 
the axial and radial perturbations were ±0.025 R. The results derived from all 
possible combinations of these axial and radial changes fill the shaded area in 
figure 14. That the above measured and computed results are roughly centered on this 
shaded region indicates that for this case the measured trajectories are fairly accu- 
rate. However, the best comparison with the measured loading occurs when the vortex 
radial location (at Y > 180") is increased (that is, the contraction is decreased) 
by 0.025 R. The identical situation was found to occur in computations of the 12" 
pitch cases; that is, the best comparison occurred when the radial vortex location 
was increased by 0.025 R over the measured value (fig. 15). For the 5" collective- 
pitch case, the situation was a little different in that a reasonable comparison of 
computation and loading data could not be made until the vortex strength was reduced 
to A ■ 0.75. In this case, the vorticity which would otherwise have been in the tip 
vortex was now included in the vortex sheet model. (For a complete description ?f 
the assumed vortex sheet model see ref. 6.) This result is consistent with the mea- 
sured vortex strength and gives the comparison sho%m in figure 16. 
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The previous comparisons have been made at low tip riach numbers. The lifting 
surface code used should be applicable to predict the spanrlee and chordvlse loading 
up to the onset of supercritical flow. Beyond this point, linear aerodynamics are 
not applicable on the blade and a more complete flow description is required. As a 
preliminary evaluation of the high-speed flow data, two-dimensional computations were 
made of the flow at the 80Z radial station. This was done using Holst's full-potential 
code (ref. 17). In order to perform this computation, an angle of attack is required. 
Since the region of supersonic flow Is localised (l.e., limited to the ins»dlate 
vicinity of the upper blade surface), it should be possible to find the angle of attack 
using the linear lifting surface code. Of course, the lifting surface code requires 
the measured vortex location and strength as mentioned previously. With the angle of 
attack obtained thereby, the Holst code produced the results shown in figure 17. This 
figure shows two computed results — an inviscid result and one with a viscous ramp- 
boundary layer model (ref. 18). It is seen that a shock-boundary layer Interaction 
model is very necessary and in this case very effective. 


5. CONCLUDING REMARKS 


The present study was intended as a benchmark to aid in the development of hov<;r 
performance codes. The goal was co obtain simultaneous measurements of blade load 
distribution and tip vortex geometry and strength using fairly standard techniques. 

In spite of some uncertainties (due mainly to wake unsteadiness), lifting surface 
computations show that the present measured loads and wake measurements are generally 
consistent with each other. 

The main conclusions from this study are: 

1. The Cook vortex measurement technique seems to be quite effective for two- 
bladed rotors. 

2. At low rotor speeds, an untwisted, untapered, double-blade-i rot-.r produces 
tip vortices which can closely resemble a classical Ranklne vortex. l^r<ci‘.pt ^or the 
lowest pitch settings, this vortex strength closely approaches the ;aaT^'.nau blade bound 
circulation. At higher tif speeds, the inner vortex structure appear'^ iocreasingly 
nonclass ical; however, the strength is unaltered. 

3. It is not possible to predict the blade-spanwlse-load di *tr ' ,.tr.-in without 

accurate vortex location and strength data. The present measur.d vorux iocation 
data were significantly different (for presently unknown reasons’ ror. the .'.l.:,sslcal 
data in the literature. However, these measurements were co obtaining 

a reasonable comparison of theor- and experiment. 

4. For the present rotor and speed range tested-. i :e op.^t.' of car^io ic flow was 

found to have no effect on the spanwise loading distr.butio*.'. ar^d tli*' vc. trajec- 
tories. The chordwlse loading is profoundly altered by the and can 

only be simulated by nonlinear aerodyiuualc technlqiv:? whl;’' 35'** .v /.lork-boundary 
layer interaction model. 
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TABLE 31.- LOCAL PRESSURE COEFFICIENT 
e “ 12* 0 - 2074 rpm M . - 0.723 
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TABLE i2.- LOCAL PRESSURE COEFPICIEHT 
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Figure 1.- The model and ex(>erimental set-up. 
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Figure 2.- Comparison of measured pressure distributions at r/8 “0.8 from each 
blade; collective pitch 9^, ■ 8“ (solid line - right blade, open symbol ■ left 
blade) • 
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Figure 4.- Measured pressure 
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TEST DATA, fig - 8“ 

• = 1 250 rpm. Cj » 0.00460 

o J2 = 2050 rpm, C j ■ 0.00461 
□ n » 2500 rpm, Cj • 0.00464 
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Figure 6.- Effect of rotor speed on blade span loading. 
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• PRESENT TEST DATA, 
dg - 12“, Mr - a22S. r/R - a96 
O GRAY'S TEST DATA, 

fig - 1 1.4*. My - a280. r/R - a966 
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Figure 7.- Comparison of present results with single blade tip loading data 
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Figure 8.- Typical wake probe data. 
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VORTEX POSITION 




Figure 


.2H 


.3 - 
.4 . 


O n = 650 rpm 
A n= 1250 rpm 

n * 2540, 2414 rpm 
□ n = 2250 rpm 

CURVE FIT OF PRESENT DATA 

KOCUREK, TANGIER DATA. Cx=0.0046 




Wake geometry measurements for various rotor speeds and comparison with 
classical data; collective pitch 0(. - 8“. 
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Figure 13. Typical velocity-time trace and 1/R curve fit for various rotor 
speeds; collective pitch 0j, - 8®, vortex age ■ 50°-65®. 
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Figure 14.- Effect of vortex position on loading computation. 
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• TEST DATA 6^ * 12°. J2R - 150 m/$ (491 fps), Cj - 0.0079 



RADIAL STATION. r/R 


Figure 15.- Compari&on of measured and computed loading. 
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• TEST DATA. 0 ^ = 5°. «R =* 150 m/s (491 fps), 



Figure lb.- Comparison of measured and computed loading. 
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Figure 17.- Comparison of measured and computed chordwise pressure distribution 



